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Abstract 



In this paper we consider a class of systems of two coupled real scalar fields 
in bidimensional spacetime, with the main motivation of studying classical or 
linear stability of soliton solutions. Firstly, we present the class of systems 
and comment on the topological profile of soliton solutions one can find from 
the first-order equations that solve the equations of motion. After doing that, 
we follow the standard approach to classical stability to introduce the main 
steps one needs to obtain the spectra of Schrodinger operators that appear in 
this class of systems. We consider a specific system, from which we illustrate 
the general calculations and present some analytical results. We also consider 
another system, more general, and we present another investigation, that 
introduces new results and offers a comparison with the former investigations. 
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I. INTRODUCTION 



Lagrangian systems described by coupled scalar fields are gaining renewed attention re- 
cently. In the case of real fields, in particular, the work presented in flj-Q has introduced a 
specific class of systems of two coupled real scalar fields. In distinction to older models, the 
class of systems introduced in the above papers is very peculiar, at least in bidimensional 
spacetime, where it presents the following general properties: firstly, the corresponding equa- 
tions of motion are solved by field configurations obeying first-order differential equations; 
secondly, the classical configurations that solve the first-order equations present minimum 
energy and are classically or linearly stable. On the other hand, the first-order differential 
equations can be seen as a dynamical system, and so we can take advantage of all the math- 
ematical tools available to dynamical systems to deal with the set of first-order equations 
and, consequently, with solutions to the equations of motion. 

For the specific issue concerning classical or linear stability, in a general way of 
investigating stability was presented. In this case, the investigation relies essentially on 
proving that the associate Schrodinger operator is positive semi-definite. This procedure of 
implementing stability investigations [^] is in distinction to the so-called standard approach 
where one usually obtains the complete set of eigenvalues of the Schrodinger operator. 

However, we know that the standard approach is important since it sets forward results 
one needs to implement quantum corrections, since in this case one has to know explicitly 
the spectrum of the corresponding Schrodinger operator [|,^. We then think that a stan- 
dard investigation of classical stability of systems of two real scalar fields belonging to the 
class of systems already introduced in is welcome, not only to present comparison with 
former investigations but also to unveil the full spectra of the corresponding Schrodinger op- 
erators. Within this context, the main motivation of the present paper is to deal with issues 
concerning finding spectra of Schrodinger operators that naturally appear in the standard 
approach to linear stability of systems belonging to the class of coupled fields introduced 
in [0-^. As we are going to show, in this class of systems the Schrodinger operators can 
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be written in terms of first-order operators, and this property eases the calculation toward 
obtaining the corresponding energy spectra. 

Evidently, to implement the standard approach to linear stability we have to deal with 
specific systems. However, to make the present investigation as general as possible, we have 
organized the present work as follows. In the next Sec. y we present the class of systems 
of two coupled real scalar fields and comment on the topological profile of the classical 



configurations. In Sec. |ITl| we investigate linear stability within the standard approach, and 
we shed new light on the issue concerning unveiling the full spectra of Schrodinger operators 
that appear in the class of systems introduced in the former section. In order to illustrate 



the general calculations, in Sec. |I^ we examine a particular system of two coupled fields. 
We end this paper in Sec. 0, where we present another investigation, in which we deal with 
a more general system, not belonging to the class of systems we shall introduce in Sec. ||. 
Each section contains some comments and conclusions. 



II. SYSTEMS OF TWO REAL SCALAR FIELDS 

Let us start with the Lagrangian density 

c = U^<pd^<p + U^xd''x-u, (1) 

where U = U{(j), x) is the potential, in general a nonlinear function of the two fields (p and x- 
Here we are using natural units, and so h = c = 1, and the metric is such that = (t, x) and 
Xa = {t, —x). The class of systems introduced in is defined by the following potential 

U{<f>,x) = lHl + ^Hl (2) 

where H = H{(f), x) is a smooth but otherwise arbitrary function of the fields and and 
H^, = dH/d(j), = dH/dx- 

The Euler-Lagrange equations, that is, the equations of motion that follow from the 
above system are given by 
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^2 



(3) 



and 



(4) 



and for static field configurations they change to 



and 



dx'^ 



fx 

dx"^ 



(5) 



H^H^^ + HyH- 



(6) 



Recall that static field configurations are configurations written in their rest reference frame. 
For static field configurations, the energy can be cast to the form E = Em + E', where 



Em = H{<Pioo), x(oo)) - i/(0(-oo), x(-oo)), 



(7) 



and 



E' 



dx 



— Ha.. 



+ 



'dx 
dx 



2n 



(8) 



As we have already learned from we impose the conditions 



dx 



dx dx ^' 



(9) 



and in this case we see that the energy gets to its lower bound Em, and the above first-order 
equations (H) solve the corresponding equations of motion (H) and (^. 

Before investigating classical stability, let us first comment on the issue concerning topo- 
logical properties of soliton solutions. From the above calculations on energy of the corre- 
sponding static fields, we see that field configurations obeying the pair of first-order equations 
present minimum energy that only depends on the difference between the two asymptotic 
behaviors of H{(j),x)- In this case, classical pairs of static field configurations having finite 



but nonzero energy must necessarily obey the topological property of connecting distinct 
minima of the corresponding potential. As one knows, this type of field configurations are 
named topological solutions 0. On the other hand, one also knows that systems of two 
coupled scalar fields may present nontopological solutions and in this case the field con- 
figurations must have the same asymptotic behavior. In the above class of systems, however, 
static field configurations obeying the first-order differential equations (P) can not have non- 
topological profile, because this would give zero energy to the pair of solutions, and zero is 
the energy value of the vacua states. This result is very interesting, because it leads us to 
the fact that the search for classical solutions of the equations of motion via the first-order 
equations (|^) evidently does not give us the full set of physical solutions, and certainly no 
nontopological solution. 

Another interesting point relies on considering the set of first-order equations as a dy- 
namical system. This procedure allows unveiling the full set of singular points, together with 
stabihty properties of each one of these points. As we can see from the potential (0), the set 
of singular points is identified with the vacuum manifold of the corresponding field theory. 
Then, once we know the vacuum manifold and the classification of each one of its points as 
stable, unstable and saddle points, we have everything one needs to deal with finding exphcit 
soliton solutions. The route toward finding explicit solutions can, for instance, follow the 
trial orbit method introduced in [0. However, since here we are dealing with first-order 
equations this trial orbit method becomes easer to be implemented than it is in the original 
work 0. 

III. CLASSICAL STABILITY 

Let us now focus attention on the issue concerning classical or linear stability. In this 
case we consider cf) = (f){x) and x = x{^) ^ P^^^ of static solutions to the above first- 
order differential equations. Here we consider 4>{x,t) = (pi^x) + rjn{x) cos{w„t) and x(x,t) = 
x{x) + ^n(a^) cos{Wnt) in order to get, from the equations of motion (|^) and (^, working up 
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to first order in the fluctuations 



This is a Schrodinger equation, and 5*2 is the Schrodinger operator, given by 



where the potential V has the form 



, V21 V22 , 

and the matrix elements can be written as 



V= , (12) 



Vu - Hj^ + H"^^ + H^H^^^ + H^H^^^, (13) 



and 



V12 - V21 - H^^H^^ + H^^H^^ + H^H^^^ + H^H^^^, (14) 



and 



V22 = H^^ + + H^H^^^ + H^H^xx- (15) 

In the above expressions a bar over H means that the corresponding quantity has to be 
calculated at the classical static values = 0(x) and x = x{^)y so Vij = Vij{x), for 

2,J = 1,2. 

On the other hand, as it was already shown in 0J^, from the first-order equations we 
can introduce the first-order operators 

+ (16) 

where v is given by 



(17) 
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Here it is not hard to check that these first-order operators are adjoint of each other, and 
that ^2 = SiS{. This is important because it proves that the Schrodinger operator 5*2 is 
positive semi-definite, and this is the result one needs to ascertain that the pair of solutions 
(f) and X is classically or linearly stable. 

To investigate the spectrum of the Schrodinger operator 5*2 we have to go further into 
this problem, and here the main difficulty concerns solving the Schrodinger equation (|10|). 
The task is not immediate since in the case of two coupled fields the fluctuations rj and ^ are 
also coupled, in general, and so the first issue we have to deal with concerns diagonalizing 
the second-order Schrodinger operator, a calculation to be done by finding the normal mode 
fluctuations. 

In the above class of systems of two coupled fields, however, we take advantage of the 
presence of the corresponding first-order operators S^, and so the task of finding the normal 
mode fluctuations is greatly simplified. This is so because here we can deal with the simpler 
task of just diagonalizing the v matrix which appears in the first-order operators. In this 
case the result allows writing 

^dx^i^Q ^_)' ^^^^ 
where the diagonal elements are given by 

v± = \{H^^ + ^xx) ± [(l/4)(^00 - H^xf + Hlx]'"- (19) 

If we use the notation ri± for the pair of normal mode fluctuations, then we have to deal 
with the following Schrodinger equations 

4±r/± = 7i;V, (20) 
where the Schrodinger operators are now given by 

St-~^V., (21) 

with the potentials 



14 = 4 + ^. (22) 

Before paying attention to specific systems, let us reason a little more on the issue 
concerning linear stability. As we can see from the above investigation, to unveil the spectra 
of we recognize that the square root that appears in v± complicates the calculation, and 
will certainly require numerical investigations. To circumvent this difficulty and perhaps 
give explicit analytical results, we should focus attention on avoiding the square root in 
v±. Here we see that the simplest case where the square root in v± desapears is when 
H^^ = 0. However, since we are dealing with systems of two coupled scalar fields, to work 
with nontrivial systems we must have H^f,^ nonzero, in order to account for interactions 
between the two fields. In this way, to get to H^^ = 0, the pair of classical solutions (f) and 
X must be very specific. However, in general the quantity H^^ does contribute, and so we 
must first deal with the quantity 

R = ^{H^^ - H^^f + Hl^, (23) 

which appears inside the square root in v±. Here we should work out a way of avoiding the 
classical pair and x to remain inside the square root. This reasoning will become clearer 
in the following, where we deal with a specific system. 

IV. AN EXAMPLE 

As a particular example, let us consider the system defined by the function 

H{<p, x) = A - a'<p^ + (24) 
In this case the potential that specifies the system has the form 

f/(0, X) = - a'f + - a')x' + \^^\' + \^^^\\ (25) 

Our main motivation to work with the above system is that it is similar to the model already 
investigated in 0, for which a standard stability investigation was already done and that 
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it presents pairs of soliton solutions [jT[ that are very similar to the pairs considered in 0. 
This motivation broadens with the fact that this specific system was shown to be useful not 
only in field theory but also in condensed matter [Q. 
In this case the first-order equations become 

^ = A(0^ - a^) + \^x\ (26) 

and 

| = .^X. (27) 

This system was already investiagted in [|l|, and some pairs of soliton solutions were pre- 
sented. In particular, a pair of solutions is 

<^i(x) = — atanh(Aax), = 0- (28) 

We recall that this pair of solutions introduces no rectrictions on the two parameters A an 
/i. Another pair of solutions is 

</'2(a^) = — atanh(/iax), (29) 
X2(a;) = 2 ^— — 1^ sech(/iax). (30) 

For this second pair of solutions the parameters A and /i are restricted to satisfy \/ ^> 1. 
Note that the limit X/^ 1 transforms the second pair of solutions into the first one. Note 
also that for the second pair of solutions the field configurations obey 

+ \l = a'. (31) 

We see that both pairs of solutions connect the points (a, 0) and (—a, 0) in the (0, x) plane, 
the first by a straight line, and the second by a elliptical line, as shown in Fig.l. 



Fig.l. The two pairs of soliton solutions. 
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As we can see, these two pairs of soliton solutions belong to the same topological sector, 
and present the same energy [0,^. Furthermore, the solutions are very similar to the pairs 
of classical configurations investigated in [p, and so it seems interesting to compare the 
present calculations with the ones there introduced. Here we recall that we already know 
that the above pairs of solutions are stable 0,^, while the pairs considered in |Q] were shown 
to be unstable, at least in the region of parameters there considered. This is a interesting 
result, since it shows that, in distinction to the class of systems defined by the function 
H = H{(f),x), older systems like the one presented in may perhaps have no classically 
stable soliton solutions, and this is true at least in some region in parameter space. 

To obtain the spectra of the corresponding Schrodinger operators we use H{(f), x) given 
by to write iJ^^ = 2A0, H^^ = ficj), and H^^ = nx- For the first pair of solutions we 
get Hfj,^ = 0, and so the equations for the fluctuations are already decoupled. This is very 
specific, and appears because of the classical value x = 0- this case we get 

f+ = — 2Aatanh(Aaa;), (32) 

and 

= — yua tanh(Aax). (33) 

For the second pair of solutions the fluctuations are coupled, and so firstly we consider 
R as given by Eg. (p3D . Here we have 

R = f^'{j^-lJ<t>l + ^^'xl (34) 
We use the orbit (^) to rewrite the above quantity as 

In this case the parameters obey A//i > 1, and so the only way of getting rid of the classical 
field from the square root is by setting to zero the coefficient of 0. Here we get the relation 
A = (3/2)/i, and this is a very interesting point in parameter space since it equals the 
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amplitude of the two classical configurations, as we can immediately see from the second 
pair of solutions; also, the value A = (3/2)/i changes the in general elliptical profile of the 



orbit (pi]) to the very particular case of a circular one. 

For the first pair of solutions, the Schrodinger operators corresponding to fluctuations 
about the and x fields are given by, respectively, 

-1) = -f- + AX^a^ - Q\^a\ech\\ax), (36) 

and 

S'o^'^'' = — — r + u^a^ — u(A + u)a^sech^(Aax). (37) 

This system was already investigated in From the results there obtained we see that, 
for A = (3/2)/i, fluctuations about the field presents the zero mode, and a bound state at 
the value w'^ = (3/4)9yU^a^, with the continuum starting at 9yU^a^. For fluctuations about 
the X field, only the zero mode is present, and the continuum starts at the value fi'^a^. 

For the second pair of solutions, the Schrodinger operators corresponding to fluctuations 
about the normal modes can be written as 

5*2^'^'' = — — + =F 4yU^a^ tanh(/iax) — 6/i^a^sech^(yuax). (38) 

dx'' 

In this case we see that these fluctuations have zero modes and no other bound states. 
Furthermore, the continua start at /i^a^, and are formed by reflecting states in the interval 
/i^a^ and and free states for energies greater than 

Before ending this paper, we notice that fluctuations about the first pair of solutions are 
described by reflectionless potentials, for which the continua start at 9/i^a^ and at yU^a^, as 
shown in Fig. 2. 





Fig.2. Potentials for the first pair of solutions. 
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On the other hand, fluctuations about the second pair of solutions are described by 
potentials that accommodate reflecting states, and these reflecting states appear in between 
the values /i^a^ and 9fi^a^, as depicted in Fig. 3. 




Fig. 3. Potentials for the second pair of solutions. 

This property in fact independs of the particular ratio between the two parameters: 
for X/fi > 1, it is not hard to show that fluctuations about the second pair of solutions 
are always described by potentials that accommodate reflecting states in between /i^a^ and 
4A^a^, which are exactly the values where start the continua of the reflectionless potentials 
that describe fluctuations about the flrst pair of solutions. 

V. A GENERAL SYSTEM 

Let us now consider another system, deflned by the potential 

U{4>, x) = \{^'- 1')' + \fx' + \\X' + \d{<i^' - 1)X'. (39) 

This is the pontential considered in and here we use the same notation, with A, / and 
d real and positive. In this case we see that 

W,O) = ^(0^-1)^ (40) 

and so this system presents the pair of solutions 

03(x) = tanh(a;/v^), Xs = 0. (41) 

In 1^ it was shown that this pair is unstable, at least in some region in parameter space. 
We recall that the investigation carryed out in was mainly concerned with stability of 
another pair of solutions, of the same form of the one given by the second pair of solutions 
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of the former example. In that investigation, it was also shown P] that when the parameters 
of the system allows for a normal mode diagonalization that leads to analytical results, the 
pair of solutions given by the above Eq. (|4T| ) is classically unstable. 

Our main interest here is simpler, and concerns investigating if there is some region 
in parameter space where the above pair of solutions is classically or linearly stable. In 
this case, the Schrodinger operators corresponding to small fluctuations rin{x) and C,n{x) 
described by 4>{x,t) = 03 + rin{x) cos{wnt) and t) = Xs + Cn{x) cos(w„t) can be cast to 
the following forms 

5f = --^ + 2-3 sech2(a;/y2), (42) 
'2) = _A^ + f-d sech2(x/y2). (43) 

The first operator presents two bond states: The zero mode and another bound state at 
w"^ = 3/2. However, for the second operator we see that a set of n = 0, 1, 2, ... bound states 
can appear, where n obeys 



n < -[Vl + 8d - 1]. (44) 

This shows that the number of bound states depends only on the parameter d, and that 
there is at least one bound state for d > 0. The energy of the bound states are given by 



f -l\Vl + 8d -l-2nl^ (45) 



To avoid instability we focus attention on the deepest bound state: Here we see that to 
ensure stability we have to impose the restriction 

1 



f >-[l + 4:d- Vl + 8d\ , (46) 

and this shows that there is room to choose d and /, keeping the corresponding pair of 
solutions stable. Furthermore, we remark that since 1 + Ad is always greater than ^/l + 8d 
for d > 0, one can not set / — > because this would introduce instability, unavoidably. 
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On the other hand, stabihty does not impose any restriction on the sign of d — f, and 
this allows going a little further on this issue since for — / > the above potential presents 
other minima. To see this explicitly, let us note that the potential (pQ]) also gives 



(47) 



Here we see that for d — f < the points {(pl = 1, % = 0) are the only possible global minima 



of the potential. However, for d — f > there are other minima, at = and 



2 d- f 



(48) 



and these minima can be local or global minima, depending on the value of the other 
parameter, A. For simplicity, let us suppose that A > (d — /)^. In this case the points 
(±1, 0) are always global minima, but when X = {d — fY there are also global minima at 

(o,±v'i/(rf-/)). 

For d — f = r > 0, and for A > we can get another pair of solutions, of the same type 
of the former one, given by 



X^i^) = tanh x 



(49) 



To investigate stability we proceed as before: Here the Schrodinger operators are 



^(4,1) 
^2 



2 

(4,2) 



We see that S'^''^^ is like 82^'''' , that is, it presents the zero mode and another bound state 



dx'^ 
d^ 
dx"^ 

.(3,1) 



+ 



d 



r — 1 



— — r sech^ 
A 



r/2 X 



+ 2r — 3r sech^ 



r/2 X 



(50) 
(51) 



at w 



(3/2)r. However, for the other Schrodinger operator the number of bound states 



are now controlled by 



n < 



1 + 



A 



(52) 



and so there is at least one bound state. From the energy of the deepest bound state we 
get, to ensure stability of the corresponding pair of solutions. 
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d 1^1 
A r ~ 4 



d / d 



(53) 



This restriction implies that A < dr. Now, if we write A = sr we obtain 



1 < ^ < - = (54) 
r d- f 

and this imphes that there are many possibihties of choosing s without destroying stabihty 
of the corresponding pair of solutions. In particular, we can choose s = 1 to write A = = 
{d — /)^, which shows that in this case the potential (|39|) presents four degenerate global 
minima: Two at x = and 0^ = 1, and two at = and = 1/r = l/{d — f). 

The above results show that when the system presents global minima at the four points 
(±1, 0) and (0, il/i/r) we can have stable solitons joining the minima (±1, 0) by a straight 
line with x = or the minima {0,±l/^/r) by another straight line with = 0. This is 
interesting, at least within the context of searching for defects inside defects, as recently 
considered in |jlO[, in the case of systems of the type considered in this section, and in 0, 



in the case of systems belonging to the class of systems introduced in Sec. |I|. But this is 
another issue, which is presently under consideration. 

ACKNOWLEDGMENTS 

We would like to thank F. Moraes for stimulating discussions, and for critically reading 
the manuscript. We also thank F. A. Brito for some interesting comments. 



15 



REFERENCES 

[1] Bazeia D, dos Santos M J and Ribeiro R F 1995 Phys. Lett. A 208 84 
[2] Bazeia D and Santos M M 1996 Phys. Lett. A 217 28 
[3] Bazeia D, Ribeiro R F and Santos M M 1996 Phys. Rev. D 54 1852 
[4] Bazeia D, Ribeiro R F and Santos M M 1996 Phys. Rev. E 54 2943 
[5] Jackiw R 1977 Rev. Mod. Phys. 49 681 

[6] Rajaraman R 1982 Solitons and Instantons (Amsterdam: North-Holland) 

[7] Rajaraman R 1979 Phys. Rev. Lett. 42 200 

[8] Boya L J and Casahorram J 1989 Phys. Rev. A 39 4298 

[9] Morse P M and Feshbach H 1953 Methods of Mathematical Physics (New York: 
McGraw-Hill) p. 1650 

[10] Morris J R 1995 Phys. Rev. D 51 697 



16 



